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Abstract 

Error estimates with optimal convergence orders are proved for a stabilized Lagrange-Galerkin scheme 
for the Navier-Stokes equations. The scheme is a combination of Lagrange-Galerkin method and Brezzi- 
Pitkaranta’s stabilization method. It maintains the advantages of both methods; (i) It is robust for convection- 
dominated problems and the system of linear equations to be solved is symmetric, (ii) Since the PI finite 
element is employed for both velocity and pressure, the number of degrees of freedom is much smaller 
than that of other typical elements for the equations, e.g., P2/P1. Therefore, the scheme is efficient espe¬ 
cially for three-dimensional problems. The theoretical convergence orders are recognized numerically by 
two- and three-dimensional computations. 


1 Introduction 

The purpose of this paper is to prove the stability and convergence of a stabilized Lagrange-Galerkin scheme 
for the Navier-Stokes equations. The scheme is a combination of a Lagrange-Galerkin (LG) method and 
Brezzi-Pitkaranta’s stabilization method 0- It has been proposed by us in USED and, to the best of 
our knowledge, it is one of the earliest works which combine the two methods, Lagrange-Galerkin and 
stabilization. Optimal error estimates are shown for both velocity and pressure. 

The LG method is a finite element method embracing the method of characteristics. The LG method has 
common advantages, robustness for convection-dominated problems and symmetry of the resulting matrix, 
which are desirable in scientific computation of fluid dynamics. Many authors have studied LG schemes 
for convection-diffusion problems ||4l[9] | 11 1[20]|2 1I I and for the Navier-Stokes, Oseen and natural convection 
problems ITl3l5llT4llT7lJl9l24l . see also the bibliography therein. The convergence analysis of LG schemes 
for the Navier-Stokes equations has been done by Pironneau Ql and improved by Siili B24| . The analysis 
has been extended to a higher-order time scheme by Boukir et al. 0 and to the projection method by 
Achdou and Guermond 0. While in these analyses they use a stable element satisfying the conventional 
inf-sup condition U3| , we extend the convergence analysis to a stabilized LG scheme. The reason to use 
the stabilized method is to reduce the number of degrees of freedom (DOF). In fact the cheapest PI element 
is employed in our scheme for both velocity and pressure, which is based on Brezzi-Pitkaranta’s pressure- 
stabilization method. Hence, the number of DOF is much smaller than that of typical stable elements, e.g., 
P2/P1. As a result, the scheme leads to a small-size symmetric resulting matrix, which can be solved by 
a powerful linear solvers for symmetric matrices, e.g., minimal residual method (MINRES) Il2l l22ll . It is, 
therefore, efficient especially in three-dimensional computation. 

In LG schemes the position at the previous time t" 1 of a particle is sought along the trajectory, which is 
governed by a system of ordinary differential equations. The position at t n ' 1 of a particle at a point at t n is 
called upwind point of the point or foot of the trajectory arriving at the point. While the system of ordinary 
differential equations is assumed to be solved exactly in 0 12411 . approximate upwind points are computed 
explicitly without assuming the exact solvability of the ordinary differential equations in |5l lT9l . Therefore, 
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we may say that the latter schemes are fully discrete. Our scheme is also fully discrete since the approximate 
upwind points are simply obtained by the Euler method. In fully discrete schemes, however, it is not obvious 
that the approximate upwind points remain in the domain, which should be proved. Such difficulty caused 
by the nonlinearity of the Navier-Stokes equations is overcome in the proof by mathematical induction, 
which has been developed in ||5ll24l . Thus, the stability and convergence with optimal error estimates 
are proved for the velocity in //'-norm and for the pressure in L 2 -norm (Theorem |T]) and for the velocity 
in /_ 2 -norm (Theorem [2j under the condition Af = Of/;''/ 4 ), where d is the dimension of the space. This 
condition is caused from the nonlinearity of the problem and it is not required for the Oseen problems ED- 
A stabilized LG scheme with an Zr-type pressure-stabilization for the Navier-Stokes equations has been 
proposed in 01, where the exact solvability of the ordinary differential equations is assumed for upwind 
points. The optimal error estimates are proved under a strong stability condition Af = 0(h 2 ) for d = 2. 

This paper is organized as follows. Our stabilized LG scheme for the Navier-Stokes equations is pre¬ 
sented in Section 0 Main results on the stability and convergence with optimal error estimates are shown in 
Section Q] and they are proved in Section[4] The theoretical convergence orders are recognized numerically 
by two- and three-dimensional computations in Section[5] Conclusions are given in Section^ In Appendix 
two lemmas used in Section 4 are proved. 

2 A stabilized Lagrange-Galerkin scheme 

We prepare function spaces and notation to be used throughout the paper. Let £2 be a bounded domain in 
R ll (d = 2,3), r = d£2 be the boundary of £2, and T be a positive constant. For an integer m > 0 and a 
real number p € [l,°o] we use the Sobolev spaces W" hp {£2), W 0 l '°°(£2), H m (£ 2)(= W m,2 {£2)), Hq(£2) and 
H~ l (£2). Foranynormed space A with norm || • ||x, we define function spaces C([0, T]\X) and //"'(0, T\X) 
consisting of A-valued functions in C([0,7’]) and //"'((), T), respectively. We use the same notation (•,•) 
to represent the Lr{£2) inner product for scalar-, vector- and matrix-valued functions. The dual pairing 
between X and the dual space X' is denoted by (•,•). The norms on W m,p (£2) d and H m (£2) d are simply 
denoted as 


II ■ ||m,p — || • II w m ’P(n) d ’ II ■ II»1 — II ■ \\H m (£2) d ( — II ' IU,2) 

and the notation || • || m is employed not only for vector-valued functions but also for scalar-valued ones. We 
also denote the norm on //“'(/ 2) d by || • ||_i. L ( 2 (£2 ) is a subspace of L 2 {£2) defined by 

L 2 (Q) = {qeL 2 (Qy, (4,1) =0}. 

We often omit [0, T], £2 and/or d if there is no confusion, e.g., C(L°°) in place of C([0,7’];L“(/2) rf ). For to 
and f[ S 1 we introduce function spaces 

Z m (/ 0 ,/t) = {v G H j (Q) d y, 7 = 0,--- ,m, ||v||zm (ro>fl) < °o} 

and Z m = Z m (Q,T), where the norm ||v|| Z m(, () ;| ) is defined by 

fm 1/2 

IMIz m (I 0 ,f 1 ) = | H ll v ll( 


We consider the Navier-Stokes problem; find (u,p) : £2 x (0, T) —> M'' x R such that 


V{2 vD(u)}+Vp = f 

in Q. x (0,r), 

(la) 

V • u = 0 

in £2 x (0,T), 

(lb) 

u = 0 

on r x (o,r), 

(lc) 

u = M° 

in £2 , at t = 0, 

(Id) 
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where u is the velocity, p is the pressure, / : £2 x (0, T) —> R f/ is a given external force, u° : £2 —> K. rf is a 
given initial velocity, v > 0 is a viscosity, D(u) is the strain-rate tensor defined by 


Dij(u) 


1 / dui duj\ 

2 \ dxj + dxi ) ’ 


i,j = l,-" 4, 


and D/Dt is the material derivative defined by 


D d 

- = -b U ’ 

Dt dt 


V. 


Letting V = Hq(£ 2) d and Q = Lq(Q), we define bilinear forms a on V x V, b on V x Q and on 

(V x Q) x (V x Q) by 

a(u,v)=2v(D(u),D(v)), b(v,q) = -(V-v,q), sf((u,p),(v,q)) =a{u,v)+b{v,p)+b(u,q), 

respectively. Then, we can write the weak formulation of (|T]k find («./?) : (0, T) —> V x Q such that for 

t G (o,r) 


(^( ? ) ,V ) +- k/ (( m ’^)W>( v ^)) = (/(O.v), V(v,?) G V x g, 


( 2 ) 


with m(0) = u°. 

Let Af be a time increment and t" = nAt for ngffU {0}. For a function g defined in f> x (0,7’) we 
denote generally by g". Let X : (0. T) —>■ S' 1 be a solution of the system of ordinary differential 

equations, 

f = “(*’')■ ® 

Then, it holds that 


£(*(,),.) = !«(*(.>,o. 

when u is smooth. Let X(■\x.t n ) be the solution of ([!} subject to an initial condition X(t n ) = x. For a 
velocity w : Q —> let X\ (w, Af) : Q —t be a mapping defined by 

Xi (w, At) ( x ) = x — w(x)At. (4) 


Since the position Aj (u" 1 .At)(x) is an approximation of X(t n l -,x,t n ) for n > 1, we can consider a first 
order approximation of the material derivative at ( x,t n ), 




u n — u n ~ 1 o X\ {u n ~ 1 , At) 
At 


(x) + 0(At), 


where the symbol o stands for the composition of functions, 


(vo w)(x) = v(w(x)), 

for v : £2 — E f/ and w : Q —> £2. X\ (w, At)(x) is called an upwind point of x with respect to the velocity w. 
The next proposition gives a sufficient condition to guarantee all upwind points are in £2. 

Proposition 1 ( ll2Tl Proposition 1]) Let w G wj '°°(£2) d be a given function, and assume 

Af||w||i j0 o < 1. 


Then, it holds that 


X l (w,At)(Q)=Q. 
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For the sake of simplicity we assume that £2 is a polygonal (d = 2) or polyhedral (d = 3) domain. Let 
= {/if} be a triangulation of £2 (= be a diameter of K £ and h = max^ e ^ /z/f be 

the maximum element size. Throughout this paper we consider a regular family of triangulations 
satisfying the inverse assumption ||8), i.e., there exists a positive constant Oo independent of h such that 

A<oo, V/f £ ■%, V/z. (5) 

h K 


We define function spaces X/,, Af/„ V/, and £>/, by 

20, = {v/ ? e c(£2) d ; v h \ K e A (/f) rf , WK e ^}, m a = {<?/, e C(fl); e A (A), V/f € ^}, 


V h = Xh n V and Q/, = M/, D Q, respectively, where Pi (K) is the space of linear functions on K £ ,5/,- Let 
Nt = /AfJ be a total number of time steps, (5o be a positive constant and (vk be the L 2 (K) d inner 

product. We define bilinear forms on H 1 (Q) x H l (£2) and ^4 on (V x H l (£2)) x (V x Z/ 1 (f2)) by 

%i{p,q) = ^o I /! l( v /L v <7k, 

KeS h 

■s^h((u,p), {v,q)) = a(u,v)+b(v,p)+b(u,q ) - % t (p,q). (6) 

The bilinear form %/, has been originally introduced in 0 in order to stabilize the pressure. 

Suppose f £ C([0,7’];L 2 (f2) d ) and zz° £ V. Let an approximate function ujj £ V* of i/° be given. Our 
stabilized LG scheme for 0 is to find {( u'f p")}„L j C V/, x Qi, such that for n = 1, • • • ,Nt 


( u 'h~ u l 1 qX i(K } 

l-Af- ’ VA J 


+ ^ h ((u n h ,p n h ),(v h ,q h )) = (f‘,v h ), 


V{vh,qh) e V/, X Q/,. 


(7) 


Remark 1 (i) By expanding u" h and p'J t in terms of a basis ofVh and Q/ v scheme (0 leads to a symmetric 
matrix of the form 

A B t \ 

B -C) ’ 

where A, B and C are sub-matrices derived from ^(zz)j,v/,) + a(w)j, V/,), biu'^q^) and %fp n h ,qh), respec¬ 
tively, and the superscript “ T ” stands for the transposition. 

(ii) The matrix is independent of time step n and regular. The regularity is derived from the fact that 
{u'liPh) = (0,0) when m} ( _1 = /" = 0 since we have 

ttIKIIo + 2v||/2(m'/0||o + E h 2 K \\W P n h \\ 2 L2 d =0 

KeSr h 

by substituting (u^—p 1 ^) £ V), x Qi , into (vh,qti) in 0. 

(Hi) There exists a unique solution (u'fpf) if Xi(u^f l ,At) maps £2 into £2. The condition is ensured if 
Af||«r 1 ||i, 0 . < c f Proposition [7] 


3 Main results 

In this section we show the main results, conditional stability and optimal error estimates of scheme 0, 
which are proved in Section [4] 

We use the following norms and a seminorm, || • ||v,, = || • ||v = II • 111- II • He* = II • lie = II • ||o- 

( m ^ 1/2 ( 


IMIrfx) = max IKIk IMI i}jx) 

' n=Q,--,Nx m ' 


A/Ik" 
11—\ 



\\ u \\p(x) = IMk m, \p\h = \ I h 2 K (Vp,Vp) K 

T ^Ke3r h 
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for m £ {1,• • ■ ,Nt} and X = L°°{L 2), l?(12) and IIHQ). D±, is the backward difference operator defined 
by 


Definition 1 (Stokes projection) For (w,r) £V x Q we define the Stokes projection ( vv/,, f,) £ V/ t x Qj, of 
( w,r ) by 

£?h((wh,rh),(vh,qh)) = £?((w,r), (vh,qh)), V(v*,? A ) S V* x Q h . (8) 

Hypothesis 1 The solution ( u,p) of (0 satisfies u £ C([0, T]\W l,00 (£2) d ) flZ 2 fl// 1 (0, T;VC\H 2 (Q) d ) and 
p £ H 1 (0,T;QnH l (£2)). 

Theorem 1 Suppose Hypothesis\J]holds. Then, there exist positive constants h o and co independent of h 
and At such that, for any pair (h,At), 

h£(0,ho], At<c 0 h d/ \ (9) 

the following hold. 

(i) Scheme 0 with u®, the first component of the Stokes projection of (m°, 0) by 0, has a unique solu¬ 
tion (; u h ,p h ) = {(u'fpl)}^ C V h x Q h . 

(ii) It holds that 


II m /iII/"“(l”) < ll M llc(z.“) +1- 

(iii) There exists a positive constant c independent ofh and At such that 


I u h ~ 


— du 

Dtou h - 


my 


\\Ph — p\\p(L 2 ) < c(At + h). 


( 10 ) 


( 11 ) 


Remark 2 Since the initial pressure p° is not given in 0, we cannot practice the Stokes projection of 
(u°,p°). That is the reason why we employ the Stokes projection of(u°, 0) and set the first component as 
u®. This choice is sufficient for the error estimates cu and also © in Theorem\2\below. 

Hypothesis 2 The Stokes problem is regular, i.e., for any g £ Lr{£l) cl the solution ( w,r ) £VxQof the 
Stokes problem, 

^{{w,r),(v,q)) = (g,v), V(v,q) £ V X Q, 
belongs to H 2 (Q) d x H l (£2) and the estimate 

IM|2 + ||ii<Ctf|M|o 


holds, where cr is a positive constant independent of g, w and r. 

Theorem 2 Suppose Hypotheses [7] and\2\lwld. Then, there exists a positive constant c independent of h 
and At such that 

\\ u h-u\\r{H) <c(At + h 2 ), (12) 

where n/ ; is the first component of the solution of 0 stated in Theorem\J\(i). 

Remark 3 Hypothesis\2\holds, e.g., if £2 is convex in R 2 , cf fT3f. 
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4 Proofs of Theorems U and |2] 


We use c, c u and C( u .p) to represent the generic positive constants independent of the discretization parame¬ 
ters h and At. c u and c lu pj are constants depending on u and (u.p), respectively. The symbol “/ (prime)” is 
sometimes put in order to distinguish between two constants, e.g., c„ and c' u . 

4.1 Preparations 

We recall some lemmas and a proposition, which are directly used in our proofs. The next lemma is derived 
from Korn’s inequality flOl . 

Lemma 1 Let £2 be a bounded domain with a Lipschitz-continuous boundary. Then, there exists a positive 
constant CC\ and the following inequalities hold. 

||£>(v)||o< ||v||t <ai||D(v)|| 0) \/veH^(Q) d . (13) 

We use inverse inequalities and interpolation properties. 

Lemma 2 ([ 8 J) There exist positive constants a, 2 i, i = 0, • • • ,4, independent ofh and the following inequal¬ 
ities hold. 


1 Qh \h 

< <*20 \\qii lo, 

G Qhi 

(14a) 

IMIo.oo 

< « 2 t/i —rf/6 ||v /j || 1 , 

Vv/, € Vh, 

(14b) 

IMIt,- 

< a 2 2^ rf/2 |Mli, 

Vv,, G V/„ 

(14c) 

l|n*v||o,» 

< ||v||o,-. 

Vv G C(£2) d , 

(14d) 

ll n />v|k=o 

< <*23 IMI 1 , 00 ) 

VvG W l ’°°(£2) d , 

(14e) 

l|n,,v-v||i 

< <*24*||v|| 2 . 

VvG H 2 (£2) d , 

(14f) 


where II/, : C(£2) d X\, is the Lagrange interpolation operator. 

Remark 4 (i) Although the inverse assumption (0 is supposed throughout the paper, it is not required for 
the estimates (fl4ab . (fl4dl) . (fl4el > and (fi4n >. The assumption that {fiTh}h\p is regular is sufficient for them, 
(ii) The inverse inequality (I14bl l is sufficient in this paper, while it is not optimal for d = 2. (iii) We note 
<*23 > 1 - 


Lemma 3 (lfl2l Lemma 3.2]) There exists a positive constant OC 30 independent ofh such that for any h 

^h((w h ,r h ),(v h ,q h )) 


inf 


(wi„rh)ev h x Q h ( VJ> m )£V h xQ Jl ( w h, r h) \Wx Q || {Vh, Qh) || Vx Q 
Remark 5 Although the conventional inf-sup condition |TS j, 

■ , b(v h ,q h ) „* 

mf SU P . >P>0, 

< ih £ Qh Vh ev h ||va||i||9a||o 


> «30- 


(15) 


does not hold true for the pair ofVand Qi„ the Pl/Pl finite element spaces, sV\ t satisfies the stability 
inequality (115b for this pair. 

Proposition2 ((§)) (i) Suppose ( w,r ) G (V H H~(£2) d ) x (QC\H ] (£2)). Then, there exists a positive con¬ 
stant 0:31 independent ofh such that for any h the Stokes projection (vt>/, , r/,) of (w, r) by © satisfies 


11 vt>/, — w 111 , ||f/,-r|| 0 , \fh~r\h < a3ih\\(w,r)\\ H 2 xH i. (16a) 

(ii) Suppose Hypothesis\2\additionally holds. Then, there exists a positive constant OC 32 independent of h 
such that for any h 

IN-HIo < «32/< 2 ||(h’,/')|| // 2 x// i. (16b) 
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We recall some results concerning the evaluation of composite functions, which are mainly due to Lemma 4.5 
in m and Lemma 1 in 0. In the next lemma a and b are any functions in W^''(Ll) d satisfying 

A?||«|| 1 , 00 , Ar11111 ;0 o < 5i, 

where 5] is a constant stated in (i) of the following lemma. We consider the mappings X\ (a. At) and 
X\ (b,At) defined in 0. 

Lemma 4 (i) There exists a constant 5] £ (0,1) such that 

J(x)> 1/2, Vx£f2, (17) 

where J is the Jacobian det(<?Aj ( a,At)/dx ). 

(ii) There exist positive constants a 4 ;, i = O, -- ,3, independent of At such that the following inequalities 
hold. 


||f-f oXi(a,At)|| 0 < <X4oA/||a|| 0i oo||f ||b 

Vf GH\Q) d , 

(18a) 

||g-g°Xi(a,At)||_i < OLiAf a Loc| g |o, 

v g eL 2 (n) d , 

(18b) 

||g oXi(b,At)-g oXi(a,At)|| 0 < OJ 42 Ar||Z? — a||oll^ll i,~, 

v g ew l '°°{n) d , 

(18c) 

||f oXi(b,At)-g oXi(a,At)|| 0 ,i < « 4 3Ar||Z? — a||ollfII 1 , 

WgeH l (Q) d . 

(18d) 


Proof 1 Since Jjj = 8jj — Atdai/dxj, (117b is obvious. It holds that for any q £ [1,°°), p £ [1,°°], p’ with 
\/p+\/p' = 1 andg £ W l ’ qp '(Q) d 

||f oXi(fc,At) -goZi(a,At)|| 0j9 < 2\\Xi(b,At) -Xi(a,Ar)|| 0 , p<? ||Vf 

from Lemma 4.5 in /0, which implies (1 1 8ab . ( 1 1 8cb and ( 118db . For the proof of ( 118bb . refer to Lemma 1 
in fffj. 


4.2 An estimate at each time step 

Let ( ui,,pi r )(t) £ Vh x Qh be the Stokes projection of (u,p)(t) by ([8]) for t £ [0,J]. Letting 

e n h = u n h -u n h , £h=p n h -Ph, r](t) = (u-u h )(t), 


we have for n > 1 


{D M e n h ,v h ) + &? h ((e n h ,£%),(v h ,qh)) = ( R n h ,v h ), V(v h ,q h ) eV h x Q,„ (19) 


where 


r',\ = L r 


i= 1 


n 

hv 


D „ _Du n u n — u n ~ 1 o X\(u n ~~ 1 , At ) 
Rhl ~ ~Dt At 

^ = ^{ T l"- 17,1-1 At)}. 


R n h2 = {m" *o Xi(u n h l ,At)-u n 1 oX\(u n 1 ,Af)|, 

R n h 4^-^{e n h ~ l -e n h - l oX iK" 1 ^)}. 


( IT9b is derived from (0. 0 and 0. We note e® = u ® — n® and set e® = p^ — p^, where (u®,p®) is the Stokes 
projection of (n°,0) by 0. 

Hereafter, let 5] be the constant in Lemma[4] 


Proposition3 (i) Let (u° ,p°) £ ( H 2 (£2) d r\V ) x (H l (Q)f\Q) be given and assume V • m° = 0. Then, there 
exists a positive constant cj independent ofh such that for any h 


vV||D( e °)|| 0 + 



< cfi. 


( 20 ) 
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( 21 ) 


(ii) Let n £ {1,-■■ ,Nt} be a fixed number and 1 £V), be known. Suppose the inequality 

Ar||«r 1 lli,-<5 1 

holds. Then, there exists a unique solution (u^,p") £ V/, x Qi, of 0. 

(Hi) Furthermore, suppose Hypothesis\I\and the inequality 

At ||m|Ic(iv 1 .° o ) — <h (22) 

hold. Let plf 1 £ Qi , be known and suppose the equation 

b(ur\q h )-V h (pl-\q h ) = 0 1 Mq h £ Qh, (23) 

holds. Then, it holds that 

Av(v||/)(e/j)||o + ^£*1*) + ^II^a/^IIo ^^ 1 (IK! 1 IIo,°°)v||£ ) (<?/'- | )||o 

+M{\Wl l \\o,o°){At\\u\\z2( tn -ij„s j +h‘-(^—\\(u,p)\\- Hl ^ij„. H 2 xH ij + l)|, (24) 

where Aj, i = 1,2, are functions defined by 

A i {Z)=c i {Z 2 + 1) 

and Ci, i = 1,2, are positive constants independent of h and At. They are defined by (ED below. 

For the proof we use the next lemma, which is proved in Annendix IA. II 
Lemma 5 Suppose Hypothesis\J]holds. Let n £ {1, • • • , Ay} be a fixed number and u n h 1 £ V), be known. 


Then, under the conditions ED and (1221) it holds that 

IKtllo < Cuv/aF||M|| z 2 (i«-i,f")) (25a) 

Iloilo <c„(|| e r 1 |lo + ^ll(«,Fr 1 ||/r 2 xffl ), (25b) 

\\R h3 \\ 0 < Q\u h || 0 , 0 ° + 1)11 { U ’P)\\H 1 (t n - l ,t n ;H 2 xH 1 )’ (25c) 

Iloilo <c|K“ 1 ||o ! oc|kr 1 Hi- (25d) 


Proof 2 (Proof of Proposition^ We prove (i). Since and are the Stokes projections of 

(m°,0) and (u°,p°) by 0, respectively, we have 

Iloilo < ikgiii = iK-sgiit < |K-« 0 II 1 + II«°-^II 1 <2«3 i*ii(m 0 , jP 0 )||^2 X// i, 

\ £ h\h = \Ph~Ph\h < \Ph-°\h + \Ph-P°\h + \P°\h < «2o(||f“-0||o+||F°-P°||o) +*||p°||l 

< (2a 2 o«3i + 1 )^II(m 0 ,F 0 )||//2 xH i, 

which imply ED ) fora = {2y/va 31 + a/<5 0 / 2 ( 2 <* 200:31 + 1 )}|| [u°,p°)\\ H i xH i. 

(ii) is obtained from (121b and Remark{J}(iii). 

We prove (iii). Substituting (D^e'j v 0) into ( Vi,,qh ) in ( fl9l ) and using the inequality (a 2 — b 2 )/2 < 
a(a — b), we have 

\\D At el\\ 2 0 + D A ,(v\\D(el)\\ 2 0 ) +b(D At e n h ,E n h ) < (26) 

1=1 

where it is noted that X\(u n ~ l ,At) in R'j^ (i = 1,2) maps Q onto £2 by (122b . From ( 123b and 0 with 
\’h = 0 £ 14 it holds that 

K u h^h) ~ ^hiPh^h) = 0, V?/, £ Q h , 
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(27) 



for k = n—\ and n. Since is the Stokes projection of(u'\p") by ©, we have 


b(u k hl q h )-%fp k lv q h ) = b(u k ,q h ) = 0, Vq h £ Qh, (28) 

for k = n — 1 and n. m and E! imply 

b(D^e n hl q h ) - %(D At e'^q h ) = 0, \/q h £ Q h , 

which leads to 


—b ( D^t e n h , el )+%i(Dkt e'l , e'l ) — 0 


(29) 


by putting qj t = — e'l £ Qi,. Adding (1291 ) to (126b and using Lemma\5\and the inequality ab < /3a 2 /2 + 
Z? 2 /(2/3) (j3 > 0), we have 

\\D^e% + D^(vMeV\\t + ^\el\i) < X>L^> 


< 


(I ft) IMII 

1=1 7 


1=1 

-n n \\2 , 1 ( 1 . \K~ l Wl 


-7T + 


V \ ft ft 


/ f At 


” l " C "i ft H“llz 2 (t'>- 1 ,«") + ^ IK m ’/ : ’)IIc(h 2 x// 1 ) + 


v|| J D(er 1 )ll5 

,«-I||2 _|_ 1 


91 No,~ 


ft 1 


ft A t 




/or any positive numbers ft (i = 1, ■ ■ • ,4), where the inequality \\e^ 1 ||o < ||e^ 1 1| i has been used. By setting 
ft = 1/8/or i = 1, • • • ,4 in (l30l) it holds that 

5A,(viiD(4)ii 2 +|| e/ :i 2 )+i||D Ai ^ii 2 <^(ii«r i ii 2 0O +i)v|iD( e r i )ii 2 

+ c («,p) jAf|M|“ 2 ( fn -l jfn ) +h 2 (\\u'} t 1 ||o j oo+ l) (^/ll( M ’^)llH 1 ((' ! ^ 1 ,f";H 2 Xi? 1 ) + l) }' 


Putting 


C] — C i( /V, C2 — C( U pp 


(31) 


we 


obtain (El). 


4.3 Proof of Theorem [lj 

The proof is performed by induction through three steps. 

Step 1 (Setting co and hf): Let c/ and A, (i = 1,2) be the constant and the functions in Proposition [3] 
respectively. Let ai, a 2 and c* be constants defined by 


at = Ai(||m|| C ( L =.) + 1), a2 = A2 (|M| C (l"o) + 1), 

c*= ■^=exp(air/2)max|a2 / 2 ||M|| Z 2,a2 / 2 (||(M,p)|| // i (// 2 xi/ i ) + 7 ’ 1/2 ) +c/j. 

We can choose sufficiently small positive constants c o and Iiq such that 

«21 {c*(c 0 /!o /12 + h l ( f d/6 ) + (a 24 + a2\)h^ d/b \\(u lP )\\ C(H 2 xHl) ^ 
a 2 2^c^co + h^ d/4 ) + (a 24 + a3i)h l 0 ^ d/4 \\(u,p)\\ c{H 2 xHl) ^ + a 2 3h‘l/ 4 \\u\\ C{w i,^ < ft, 


< 1, 


CO 


since all the powers of /zq are positive. 


(32a) 

(32b) 
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Step 2 (Induction): For n £ {0, • • • . Nt } we set property P(n), 


( a ) v ll^ ) ( e /!)llo + "2"l e A \l + 2 \\ D ^ e h\\\(L 2 ) 


POO : 


2 2 
<exp(ai«Af) v|| D{e®) 


p 

o + ^\ e h\h + a2 ^''\\ u \\~z^(o,t n ) + / r (ll( M ’P)llr/i( 0 ,f" ; //2 x //i) + «Af) j 


( b ) IKIk- < ||m||c(l”°) +1) 

.(c) At\\u n h \\i^<8i, 


where ||DAre*||; 2 f/ 2 j vanishes for n = 0. P(n)-(a) can be rewritten as 

n ( " \ 
x„ + A t'Y_ i yi <exp(ainAr)(xo + Ar^I?,'j, 


(33) 


where 


x„ = v||D(e/|)||o+ ■ 


K\ 1 , 


yt = -||£>Af4llo) 


bi = «2i Af ||M 


I z 2 (f- 


/) 


+ /r 


vAf ^ ( M,p )ll 


if 1 


+ 1 


We firstly prove the general step in the induction. Supposing that P (n — 1) holds true for an integer 
n £ {1, • • • ,Nt}, we prove that P (n) also does. Since P (n — l)-(c) is nothing but (fTil) . there exists a unique 
solution (u1 v p'il) £ Vh x Qu of equation © from Proposition®(ii). We prove P(n)-(a). (l22l) holds from the 
estimate. 


A f ll“llc(W l ’“’) < CQhJ ||w||c(Wl’°°) < CO a 23hJ llwllcfiv 1 '”) — ^ 1 ) 

from condition @, RemarklH-(iii) and (132bl) . (123b is obtained from © for n> 2 and from the definition of 
i.e., the Stokes projection of (m°, 0) by ®, for n = 1. Hence (f24l i holds from Proposition©(iii). 
Since the inequalities ||o,oo) < at (i= 1,2) hold from P(« — l)-(b), (f24l) implies 


D A ,x n +y n < aix n -i+bn, 


which leads to 


x„ + A ty n < exp(a 1 Ar)(x„_i + A tb„) 
by 1 < 1 + cqAf < exp(oiAr). From P(n — l)-(a), i.e., 

n— 1 . n— 1 . 

x„-\ + Ar ^ yt < exp{fli(n — l)Ar} (xo + At Y bi ), 
1=1 1=1 

it holds that 

n n— 1 

x„+At Y^yt < exp(fliAr)(x„_i +Atb n )+At Y yt ( b y C©) 

i= 1 i= t 

n~ t 

< exp(fl!Ar)(x„_i +Ar Yyi + ^tb,,) 

v i=t ' 

< exp(fliAf) expjfli (n — l)Ar} ^xo + At Y b^j + A tb h 
>tp(fli«Af) ^xq + At Y bi 


< exr 


i—1 


(by®) 


(34) 


(35) 


which is (l33l >. i.e., P(n)-(a). 

For the proofs of P(u)-(b) and (c) we prepare the estimate of ||ejj||i. From P(n)-(a) and (l20l > it holds that 


v||£ , (e£)||o+ yl e £l/l + 2 H' DAfe/l ll|(i. 2 ) — ex p(aiF) cjh 2 + a 2 |Ar|| m||~ 2 +h (\\{u,p)\\hI(h 2 xh 1 ) + 7 ’)} 
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(36) 


< expffliT) 


Cl 2 At II u 


& +h 2 \ 


a 2 (11 i u > P)\\h 1 (H 2 xH 1 ) 


r)+cj}] < {c 3 (At+h)} 2 


where 

c 3 =exp(a 1 r/2)max|fl2 / 2 H«llz2. a 2 2 (W( u ’ p)\\hHh 2 xh') + 7 ’ 1/2 ) + c/}- 

( l36t implies 

\\eUi<cc l \\D(el)\\ 0 <-^c 3 (At + h)=c*(At + h). (37) 

We prove P(n)-(b) and (c). Let LI/, be the Lagrange interpolation operator stated in Lemma [2] It holds 

that 


Kilo,- < K-n A « ,, |lo,- + l|nK||o,~ < a 2l h- d / 6 \\u n h -n h u n \u + ||nK||o,- 

< « 2 i/z- rf/6 (||^ - fiJHi + ||«j - «"|l 1 + IK - n*«li) + ||n /1 M n || 0 ,<» 

< a 2 ih~ d/6 {c*{At + h) + a 3 ih\\(u’ 1 ,p n )\\ H 2 xH i + a24h\\u n \\2} + \\u”\\o,o* (by d37]») 

< a 2 l {c^c 0 hQ /l2 + h^ cl/6 ) + {a 24 + a 3 l )h^ d/ 6 \\(u 1 p)\\ C[H i xH i ) } + \\u\\ c(L ^ (by ®) 

< 1 + IMIc(l~) (by 023), 

At\\ul\\ ha o<coh d/4 (\\u n h -U h u n \\ l: o.+ \\n h u n \\ l:a o)<coh d/4 (a 2 2 h~ d/2 \\u n h -U h u n \\i + \\U h u n \\i^) 

< co{a22h~ d/ 4 {\\u n h - u n h \\i + K~ M li + \\u n -n h u n \\i) + h d/ 4 \\n h u n \\ 1 ^} 

< co[a 2 2 h~ d/4 {c*(At + h) + a 3 ih\\(u n ,p n )\\ H 2 xH i + a 2 Ah\\u n \\ 2 } + a 23 h d ^\\u n \\x^\ 

< co[cc 22 h~ d / 4 {c*(coh d / 4 + h) + (a 2 4 + a 3 i)h\\(u n ,p n )\\ H 2 xH i} + a 2 3 h d / 4 \\u n \\i } oo\ 

< co[a 22 {c t (c 0 + h l 0 ^ d/4 ) + (a 2 4 + a 3 l )hl^‘' / 4 \\(u,p)\\ c(H i xHl) } + a 23 hQ / 4 \\u\\ c{w ^ ) \ 

< 5j (by (l32bll ). 


Therefore, P(n) holds true. 

The proof of P(0) is easier than that of the general step. P(0)-(a) obviously holds with equality. P(0)-(b) 
and (c) are obtained as follows. 

Kilo,- < K - nKlk- + ||n*M°||o,- < a n h- d i\\\ul - «°||! +1|«° - nK|| 0 + HnKlk- 

< a 2 i ( 0 : 31 -(-ct24)^ 1 d/6 \\u °\\2 + ||m°||o,» < 1 + ||m||c(l“) (by ( 132a| >). 

Af||^||l^<C0/t d/4 (||MK n / ! M 0 || 1 ,~ + l|nK|| 1 ,oo)<C0/t d/4 (a22^ d/2 ||MK n ^°lll + ll n ^ 0 |ll,~) 
<co{a 22 /t“ d/4 (||^-M 0 ||i + ||M°-nK||i)+/t d/4 ||nK||i,4 

< c 0 {a 22 (a 3l + a 2 4 )h l ~ d/ 4 \\u °\\ 2 + a 23 h d/ 4 \\u°\\ hoo } < 81 (by ((320). 

Thus, the induction is completed. 

Step 3: Finally we derive the results (i), (ii) and (iii) of the theorem. The induction completed in the previous 
step implies that P(Nt ) holds true. Hence we have (i) and (ii). The first inequality of (fill in (iii) is obtained 
from d37l and the estimate 


ll M * — “llrp? 1 ) — Kllrpd) + IMIrpr 1 ) — Kllrpf 1 ) +cc 3 ih\\ ( u ,p)\\c(h 2 xh >)• 

Combining the estimate 


- „ (j Isl 

DfaUfo - 


du n 

dt 


— „ du n 

DjAjU -— 

dt 


<||D A K,||o + KA,t7lo- 

a 3 \h„, , „ I At 




d 2 i, 


dt 2 


L 2 (t"- l ,t n \L 2 ) 
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with (l36l) . we get the second inequality of ( ITTI) . Here, for the estimates of the last two terms, we have used 
the equalities 


(D At ri")(x) = j -^-(x,t n ~ l +sAt)ds, (o Ar u n - (x) 

We prove the third inequality of (ITTb . It holds that 

1 ^ h ((e n h ,e"),(v h ,q h )) _ 


|<II(«)IIvxq< 


0:30 (^)eU \\( v h,qh)\\vxQ 



+ sAt)ds. 


1 su (R’i„Vh) ~ (Dit^Vh) 

a 30 {vhAh)eV h xQ h \\(Vh,qh)\\vxQ 


< c («,p){v^ll M llz 2 (f"- 1 ,f") ^ b(y-j=\\{u,p)\\ H l^ n -l t n. H 2 xH l^ + l) + II e 'l, 'ill + HA^Hoj (38) 

for n = l,--- ,Nt- Here we have used Lemmas [3 and 0 the inequality ikr'iio < ikr'iii and ([Hi. We 
obtain the result by combining (l38l >. (l36b and the estimate 

\\Ph-p\\l2(I?) < l|£/ ! |lz2( Z 2) + ||p/ ! -.p||,2( L 2) < 1181,11,2^2) +Vfa 3 lh\\(u 7 p)\\ c{H 2 xH ly 


4.4 Proof of Theorem |2| 

We use the next lemma, which is proved in Appendix lA.21 

Lemma 6 Suppose Hypotheses [7] and [2] hold. Let n £ {l,--- ,Nt} be a fixed number and u'^ 1 £ V/, be 
known. Then, under the conditions (ED and <E3 it holds that 

Mo <c u {\\er%W\\{u, P y-'\\ H ^), (39a) 

Il^fc3 II VjJ — c u(\\{u,p) n * \\h 2 xH 1 II e 'h ' llo + —?= II ( u iP) II// 1 (/"- 1 , t n \H 2 xH l ) + ^ II { u iPY 'll H 2 xH\ 


IK 4 llv'<£«(l+^ /6 |l4“'lli)(lkr'llo + /i 2 ||(«,Pr'llrr2 x// i). 


(39b) 

(39c) 


Proof 3 (Proof of Theorem^ Since we have ||e/,||/«. (i/i) < c*(At + h) < c*(co + h l 0 d ^ A )h d / A from (ED and 
(fD, (139cb implies 

iK4iiv'<^(ikr'iio+/i 2 ii(«,pr'ii H 2 Xi ,i). (40) 

Substituting {e^—e’f) into ( 17 , ,< 57 ,) in (1191 ) and using Lemma\T\ (I25ab . d39ab . ( I39bb . (140b and the inequality 
ab < pa 2 /2 + b 2 /( 2 p) (/3 > 0 ), we have 


^(^IWIlO+SlWIli + ^l^t^) 

\ Z, / CC] j = j 


ic "(ft + 

h 4 

ftsA t 


1 ll(“’I : ’)llc(i/ 2 xff 1 ) c;\ 


+ 


f)lkr'llo+(LA)ll4ll?+4[^ll«llV 1 / 


ft ft 

IK M !/ 7 )ll|l(;"-l, f »;ff2 Xiff l ) +/l 4 {(^ + ^)||(M,ft||c( ff 2 ><if l) + ^ £ ll( M ’f’)llc ( (//2 xi/ l)} 


for any J3; > 0 (i = 1, • • • ,4), where the inequality ||e^||o < lk/|||l bas been employed. Hence, it holds that 

°A» (2 H e *Ho) + ^2 Ikftlll ^ c (u,p)\\ e 'h ' 110 + c («,p) ( 11 “ 11 z 2 1 + ^■||( M !f’)llffi(r«-i i r";// 2 x// i) +h 


by setting /3, = v/(4tt 2 ) (i = l,--- ,4). From the discrete Gronwall’s inequality there exists a positive 
constant C 4 independent of h and At such that 


ll e /ill/“(L 2 ) < c 4 (lk/!llo + Ar + h 2 ). 
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(u 0 ,p°)\\ H 2 xH l, 


Using (1 1 6hb . we have 

\\e 0 h \\o<H-A\o + h 0 -a 0 h \\o<2a 32 h 2 \ 

\\uh-u\\ nL2) < ||e A || r(L 2 ) + ||T7|| r(i 2) < ||e/,||/~(i, 2 ) + U 32 h 2 \\(u,p)\\ C ( H 2 xH iy 
Combining these three inequalities together, we get ©. 


5 Numerical results 


In this section two- and three-dimensional test problems are computed by scheme 0 in order to recognize 
the theoretical convergence orders numerically. 

For the computation of the integral 

f u n h ~ l oXi(m^ -1 ,At)(x)Vh{x) dx 
J K 

appearing in scheme 0 we employ quadrature formulae 123ft of degree five for d = 2 (seven points) and 
3 (fifteen points). The results obtained in Theorems Q] and [2] hold for any fixed bo. Here we set 8 q = 1. The 
system of linear equations is solved by MINRES 121 1221 . 

Example 1 In problem 0 we set Cl = (0, l) d , T = 1 and four values of v, 

v = 10“*, k= 4. 


The functions f and u° are given so that the exact solution is as follows: 
ford = 2 

u(x,t) = VW,), p(x,t) = sin{7r(xi+2x2 + 1)}, 

V dX2 dx i / 

\/3 

yr(x,f) = 7— sin 2 (7rxi)sin 2 (7TX9)sin{7r(xi + X9 +/)}, 

2k 


for d = 3 


u(x,t) = rot x F(x,l), p{x,t) = sin{7r(xi +2x2 +X3 + /)}, 

8-\/3 2 2 

T* i(x,f) = —— 1 sin(^xi)sin“(7TX9)sin“(^X3)sin{^(x2+X3+f)}, 
21k 

8v/3 9 2 

T , 9(x,l) = ——sin“(^xi)sin(7TX2)sin (^X3)sin{7r(x3+xi+l)}, 
21k 

Sy/3 t o 

T , 3 (x,t) = -sin“(^X! )sin 1^X2) sin (71x3) sin {n:(xi +X9 +f)}. 

21 k 


These solutions are normalized so that ||m||c(L”) = ll/ 7 llc(z.°°) = 1- 

Let N be the division number of each side of the domain. We set N = 64,128,256 and 512 for d = 2 
and N = 64 and 128 for d = 3, and (re)define h = 1 /N. Local meshes are shown in FigureQ]for d = 2 (left, 
N = 64, in [0.9, l] 2 ) and 3 (right, N = 64, in [0.9, l] 3 ). Setting A t = y\h and 72 h 2 (71 =4, 72 = 256), we 
solve Example[T|by scheme 0 with u° h , the first component of the Stokes projection of (u°, 0) by ®. The 
two relations between At and li, i.e.. At = 7 h and 72 h 2 , are employed in order to recognize the convergence 
orders of ([111 and (IT2l) . respectively and we have (At =)yih = y^lr for h = 1 /64. For the solution (uh,Ph) 
of scheme 0 we define the relative errors Erl and Er2 by 


Erl = 


\\u h -n h u\\i2 (H i) + ||p/,-n/,p|| /2(L 2 ) 

ll n A M lb 2 (ff i ) + lirvil^) 


\\u h -n h u\\ nL2) 

ll n /+llr(L 2 ) 


where for the pressure we have used the same symbol II/, as its scalar version, i.e., II/, : C( fl) —> M/,. 
Figure[2]shows the graphs of Erl versus h lord = 2 and 3 (left. At = J\h) and Er2 versus h for d = 2 (right. 
At = 7>/z 2 ) in logarithmic scale, where the symbols are summarized in TableQ] The values of Erl, Er2 and 
the slopes are presented in Table [2] We can see that Erl is almost of first order in h for both d = 2 and 3 
and that Er2 is almost of second order in h. These results are consistent with Theorems [T| and [2] 
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Figure 1: Local meshes for d = 2 (left, N = 64, in [0.9, l] 2 ) and for d = 3 (right, N = 64, in [0.9, l] 3 ). 


6 Conclusions 

A combined finite element scheme with a Lagrange-Galerkin method and Brezzi-Pitkaranta’s stabilization 
method for the Navier-Stokes equations proposed in [T5JHS1 has been theoretically analyzed. Convergence 
with the optimal error estimates of 0(At + h) for the velocity in //'-norm and the pressure in L 2 -norm 
(TheoremQ]) and of Ot At + h 2 ) for the velocity in L 2 -norm (Theorem^ have been proved. The scheme has 
the advantages of both method, robustness for convection-dominated problems, symmetry of the resulting 
matrix and the small number of DOF. We note that it is a fully discrete stabilized LG scheme in the sense 
that the exact solvability of ordinary differential equations describing the particle path is not required. In 
order to provide the initial approximate velocity we have introduced a stabilized Stokes projection, which 
works well in the analysis without any loss of convergence order. The theoretical convergence orders have 
been recognized numerically by two- and three-dimensional computations in Example Q] It is not difficult 
to consider a fully discrete stabilized LG scheme of second order in time due to the idea of BUD, and its 
convergence with the optimal error estimates will be proved by extending the argument of this paper. 
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Table 1: Symbols used in Figure[2] 




V 


d 

io-' 

\(r z 

10~ J 

10~ 4 

2 

O 

□ 

A 

V 

3 

• 

■ ▲ 

▼ 
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Figure 2: £>1 vs. h for d = 2 and 3 (left. At = 71 / 2 , 71 = 4) and Erl vs. h for d = 2 (right. At = 7 >/i 2 , 72 
256). 


Table 2: Values of Erl, Erl and slopes of the graphs in Figure [2] 




N 


Erl 


Erl 


<N 

II 

^3 

slope 

d = 3 

slope 

<N 

II 

^3 

slope 

V = 

KT 1 : 

64 

7.24 x 10 ~ 2 

— 

6.37 x 10 ^ 2 

— 

1.03 x 10 1 

— 



128 

3.85 x 10 ~ 2 

0.91 

3.25 x 10 ~ 2 

0.97 

2.96 x 10 ~ 2 

1.80 



256 

1.99 x 10 ~ 2 

0.95 

— 

— 

7.71 x 10 ~ 3 

1.94 



512 

1.01 x 10“ 2 

0.97 

— 

— 

1.96 x 10 ~ 3 

1.97 

V = 

l(r z : 

64 

1.70 x 10 " 1 

— 

2 . 10 x 10 1 

— 

2.74 x 10 ' 1 

— 



128 

9.51 x 10 ^ 2 

0.84 

1.10 x 10 _1 

0.94 

8.66 x 10~ 2 

1.66 



256 

5.13 x 10 ' 2 

0.89 

— 

— 

2.35 x 10 ~ 2 

1.88 



512 

2.68 x 10“ 2 

0.93 

— 

— 

6.09 x 10 ~ 3 

1.95 

V = 

10 ^: 

64 

2.14 x 10 " 1 

— 

3.78 x lO ^ 1 

— 

3.41 x 10 ' 1 

— 



128 

1.21 x 10" 1 

0.82 

2.02 x lO ^ 1 

0.90 

1.10 x lO ^ 1 

1.63 



256 

6.63 x 10 ' 2 

0.87 

— 

— 

3.03 x 10 ~ 2 

1.86 



512 

3.51 x 10 “ 2 

0.92 

— 

— 

7.88 x 10 ~ 3 

1.95 

V = 

kFT 

64 

2.39 x 10 ' 

— 

4.45 x lO ^ 1 

— 

3.50 x 10 ' 1 

— 



128 

1.35 x 10 ' 1 

0.83 

2.35 x lO ^ 1 

0.92 

1.13 x 10 ' 1 

1.63 



256 

7.34 x 10 “ 2 

0.88 

— 

— 

3.13 x 10 ~ 2 

1.85 



512 

3.88 x 10 “ 2 

0.92 

— 

— 

8.14 x 10 ~ 3 

1.94 
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Appendix 

A.l Proof of Lemma \ 5 \ 

Let t(s) = f" _1 + sAt (s£ [0,1]). We prove (125al >. Let y(x,s) = x— (1 — s)u n ~ l (x)At. It holds that 


*«(*) = {( 


u(y(x,s),t(s)) ^ 
d 


= At 


j^—+ m" *(x) • vjnj(x,t") + {((w" — u” ! )(x) • V)m"}(x) — J +u n 1 (x) • V^m| (y{x,s),t{s))ds 

MM ^+w" ! (x) • u}(y(x,si),t(si))dsi+AtJ |^^(x,t(i)) • vjw"|(x)ds 

1 S '{(.§-t +Un ~ l{x) - V ) U }iy( x ^i)^i))dsi+At J {(^(x,t(s)) • 


= At 

=Ku W+^ 12 W- 


Each term 7?^ is estimated as 


||**iillo<Ar J si {(^+ M " 1 (')‘ v ) M }(y(- + i)+(st)) 0 ^i <c u VAt\\u\\ Z 2 ^ t n-i r) , (A.la) 
Iloilo <c„Aty o —(.,t(s)) ds<c u y/At 


du 

dt 


L 2 (t n ~ l ,t";Lr ) 1 


(A. lb) 


where for the last inequality of (1A. 1 al l we have changed the variable from x to y and used the evaluation 
det(<9y(x,si)/<9x) >1/2 (Vsi £ [0,1]) fromLemma[4]-(i). From (IA. 1 b we get (125al) . 

(I25bb is obtained as 


Iloilo < «42||^r 1 -«”“ 1 |lo||«"- 1 || 1 ,oo < C^II^-^ll^dlTT^-^lo-Hkr'llo) 

<a* 2 ||«"- 1 || ll -(a 3 i*||(« ) p)"- 1 || /f 2 X Hi + ||« 2 - 1 ||o). 

We prove ( 125cb . Lety(x,s) =x— (1 — s)u"~ 1 (x)At. Since it holds that 

i r , , . . ,i> r l r td 


(A.l) 


Rh = ~ 


^ n(y(-,s),t(s)) s=0 = J 0 {(^7 +m /< H')■ } (y('i s )’t( s ))d s i 


we have 


ll^allo </ ) 1 |{(|: + «r 1 (-)-v)t7}(y(- a), fW) 

j ( 0 +IK“ 1 ||o,oo||VT 7 (y(-,s),f(s))|| 0 )c/s 

Vr; (-,/(s)) || 0 }js (by LemmaQjL)) 

L 2 (t"-',i":L 2 ) + ^ L ‘ h ^ 0 lOO lllli 2 (f"- 1 ,*";i 2 )) 


< 


< v 2 


^(■.>M)|| 0 +ii«r 1 i 


< 


< 


2 / d~q 
At V dt 


« 3 lMIK *110,00+1)11 (u, p)\\ H l( t n-l jn.jp xH l), 


which implies (I25cb . 
(I25db is obtained as 


WR’UWo = ~fa\\ e h l - e h * oZ lK *> A 0||o - a 4o|| 


,11 


o^ikr'iit. 
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A.2 Proof of Lemma [6] 


( 139al ) is obtained by combining (1 16bb with ( 1 A ,2b - For ( 139bl ) we divide R'j i3 into three terms, 

/ ? «3=5 Af r7» + l{T7' l - 1 - T 7 "- 1 oZ 1 ( M n - 1 ,At)} + l{r 7 ”- 1 oX 1 ( M ' ! - 1 ,AO-r 7 ' l - 1 oZi( M ”- 1 ,At)} 

— J? n I J? n 

= K h3 1 “I" K h32 + K-h33 • 

It holds that, in virtue of (1 1 6bb . 


||^3lllv'<l|OA,t?"||0<^ 


drj n 


dt 


L 2 (t n - l ,t";L 2 ' 




WR’wWvl < a 41 ||M n_1 || 1 ,oo||TI n_1 ||o < «4l|| M "- 1 || 1 ,oo a 32 h 2 Uu, P r~% 2xHl , 

\\Rm\\v; = su p |4i( j i"" loIi K" 1 ^)- j i“" loIi («"" 1 . A( )> v j) 

h v h eV h \\Vh\\l At \ ) 


(A.3a) 
(A.3b) 


< SU P TrArX7ll T? " 1 qX iK 1 ,At)-T7" 1 oI[(/ 1 ,A?)11„ j||v/j|| 0 ,oo 
v h ev h IIVftHi At 

< «43||m* _ 1 — M" _1 || 0 ||T7"“ 1 || 1 «21^ _rf/6 (A.3c) 

< « 21 «43/*- d/6 ||r 7 "- 1 ||i(|kr 1 ||o+||r7”- 1 ||o) 

<a 2 ia 43 a 3 2h 1 - d/ %u,p) n - 1 \\ H 2 xHl (\\e n h - l \\o + a 32 h 2 \\{u,p) n - 1 \\ H 2 xHl ) 

<c||(«,p)"- 1 || fl 2 xfl i(ll^" 1 ||o + A 2 ||(«,p)"- 1 ||fl 2 xfl i). (A.3d) 

From (!A.3aL (lA.3bl) and (IA.3db we obtain (I39bb . 

For (139cl) we use the estimate of R"^. R" i4 is obtained by replacing 1 with —e n h 1 in R'f l4 ~, +R \ 33 . 
Hence, from (lA.3bb and (lA.3cb we have 


||^ 4 ||V^ < « 4 1 || M "- 1 || l , oo | k"- 1 ||o + « 2 I « 43 A -^ /6 ||^- 1 || 1 || M ”- 1 - M "- 1 || 0 

< «4 I i| M "- i iii,ooik”- 1 ii 0 +« 21 «43A- rf/6 ii^- 1 ii 1 (ikr 1 iio+«32A 2 ii(«,^r- l ii^x//0 

<c„(i+/t- d / 6 i| e r i ii 1 )(i| e r i iio+/t 2 ii( M ,/ J r i ii //2xfll ), 


which implies ( I39cb . 
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